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Abstract 

We study the spectral properties of certain non-self-adjoint matrices associated with 
large directed graphs. Asymptotically the eigenvalues converge to certain curves, apart 
from a finite number that have limits not on these curves. 

1 Introduction 

In an earlier paper we studied the spectral properties of small, possibly random, perturbations of 
annxn Jordan matrix, and proved that most, but not necessarily all, of the eigenvalues lie close 
to a certain circle as « ^ oo. In this paper we consider a similar question when the perturbation 
is a sparse matrix, whose entries need not be small. 

It turns out that this problem is naturally associated with a type of directed graph that is 
relevant to a variety of problems involving unidirectional flows along one-dimensional channels 
between several junctions. Since the entity flowing may be a fluid (e.g. blood flow through the 
network of veins and arteries), traffic or electronic data, the ideas here have potentially wide 
applicability. The results in this paper also provide asymptotic spectral information for a special 
class of large directed graphs. 

The spectral theory of a directed, acyclic graph is not particularly interesting. By labelling 
the vertices appropriately one may make the adjacency matrix of the graph upper triangular, so 
its eigenvalues coincide with its diagonal elements. This correctly suggests that the interesting 
spectral properties of directed graphs are highly dependent on the the structure of its cycles. 

Let 5 be a finite set and A, y real or complex numbers indexed by /, j G S. We make S into 
a directed graph {S, —^) by putting / —>■ j if A, y 7^ 0. We will see that S can be partitioned into 
disjoint 'junctions' Ji,J2,...,Jk linked by one-dimensional 'channels' Ci,C2, ...,C/j. A natural 
way of constructing the graph is to draw the junctions as disjoint localized regions, and then 
join them by directed channels. However, one may also represent the chamiels by a sequence 
of horizontal lines, possibly of varying lengths, all directed from left to right. The junctions 
collectively may then be regarded as quasi -periodic boundary conditions that join the right ends 
back to the left. In many of the applications the entries of A are non-negative and the spectrum 
of A may also be investigated by using the Perron-Frobenius theory, IH. The theorems in this 
paper do not depend on this assumption, which is not appropriate in some contexts. 
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Given the above assumptions, we construct a new graph (5^"), — >) for every natural number 

(n) 

n as follows. We leave the junctions unaltered and replace each channel C,- by a new channel C- 
with the same endpoints but with #(c|"^) = « x #(C,). Under suitable assumptions the matrix A 
induces an associated matrix A'-"' whose coefficients are parametrized by pairs of points in S^"\ 
Our goal is to investigate the asymptotic behaviour of Spec(A(")) as « ^ oo. 

The continuous analogue of our problem replaces each channel by a continuous bounded 
interval within which the relevant operator is aD + bl, where D is first order differentiation and 
the constants a, b vary from one interval to another. Instead of junctions one has to specify 
suitable boundary conditions connecting certain groups of end points. In contrast to 13, we 
do not restrict attention to self- or skew-adjoint problems. We finally mention that the wave 
equation on an undirected graph can be modelled in these terms by replacing each unoriented 
edge by a pair of oriented edges with the same ends and associating the operators itD with the 
two edges; see Q. 

Every constant in this paper is independent of the asymptotic parameter n. 



2 Zeros of some analytic functions 

The material in this section is needed for the proofs of many of the later theorems, but does not 
mention graphs explicitly. The results that we attain are to some extent analogous to calculations 
in |[I]|l|5l|71l3 but the details are different. 

We will often need to count the number of zeros that an analytic function has in a particular 
region. The following elementary result enables us to use a change of variable to simplify this 
task. 

Lemma 1. Let U and V be open subsets ofC and g :U conformal. Iff :V is analytic 
then Fog has a zero of order mat uq €zU if and only ifF has a zero of order m at g{uo). 

In the proofs to follow we shall make repeated use of a particular type of contour which is 
perhaps best described as being the boundary of a sector of an annulus centred at the origin. The 
contour Cr,i?.e,.e2 is completely specified by four parameters r,R, 61,62 where r < R correspond 
to the two radii of the annulus and 61 < 62 correspond to the angular sweep of the sector. More 
precisely, the contour Cr,R,ei,e2 = {7i > Y2, 73; 74} is the concatenation of the following four curves 

Yi{t)=te'^' r<t<R 

Y2{t)=Re" 6i<t<62 

Y-i{t) = {r + R-t)e'^'- r<t<R 

74(0 =re'(^'+''2-') 6i<t<62 

We note that Cr.R^e^,e2 encloses the open region {z-r<\z\< R, 61 < argz < 02}- 

Lemmas [2] and [3]provide examples of the more general analysis that starts with Theorem |2l 

Lemma 2. Let p„ denote the polynomial 

Pniz) = {z- aYiz - by + a{z -ay + piz- b)" + 7 (1) 
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where a b. If a, [5 and y are all non-zero then the roots ofp„ converge as n^°oto the union 
of the two circles \z — a\ = I and \z — b\ = 1. 

Proof. Given e > put 

Se = {z:\z- a\ > l + e}n{z: \z-b\ < 1-e}. 



If z G Sg then 



Therefore 



.a = {z-bf + ^('^] + 



(z — a)" \z — aj (z — a)" 

Pn{z) 



a 



\a\ 
- ~T 



{z-aY 

for all large enough n, uniformly with respect to z G Sg. This implies that Pn{z) / 0. By carrying 
out three other similar calculations we find that every root of pn Ues in the e-neighbourhood A^e 
of the union of the two circles provided n is large enough. 

We now show that if zo lies in the union of the two circles then there exists a sequence Zn of 
complex numbers such that Zn zo and Pn{zn) = for all n. There are essentially two cases to 
consider, depending on whether zo hes on one or both of the circles |z — a| = 1 and \z — b\ = \. 

Firstly, suppose zo lies on only one of the circles. For the sake of definiteness we suppose that 
|zo — a| = 1 and |zo — ^| < 1- Let ?7 be a small open neighbourhood of zo so that | (z — a) (z — ^) | < 
c < 1 and |z - ^1 < c < 1 for all z in U. This implies that ^„(z) = (z - a)"(z -b)" + l5{z- b)" is 
uniformly exponentially small for z in i7. This suggests that inside U the roots of 

p„{z) = a{z-af + Y+qn{z) 

should be close to the roots of /„(z) = o;(z — a)" + 7. 

Suppose Zo = +a. Let w„ = \ Y/oc\^^"e"^" +a be a root of /„(z) which is closest to zo- 
Since \Y/a\^^" 1 and 0„ -> 0o we see that w„ zo- Let C„ be the contour Cr„,ff„,0i„,e2„ + « 
where r„ = (^|7/a|)i/", /?„ = (||7/o:|)i/«, 0i „ = - n/2n and 02,« = <l>n + ^/2n. The region 
enclosed by C„ contains w„. 

Let M = I7I/2. We show that z G C„ implies |/«(z)| > M, where we note the bound is 
independent of n. Since C„ = {71,72,73,74} there are four cases to consider. If z G 7i then 
z = pe'^' " + a where r„< p <Rn and so 

|/„(z)| = |apVW"-'^/2) + 7| 

= |ap"e'("^"-^/2) + 7- {air/ale'"^" + r)\ 
= l«lk'(-P")-|7/«ll 

> in 

where the second equaUty uses a 1 7/0: 1^'"*^" + y = = 0. The case z G 73 is similar. If z G 
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then z = Rne'^ + a for some 6 and so 



> - 



\Mz)\ = a-\r/a\e"" + y 

l\7\-\7\ 
= \\7\ 

The case z G 74 is similar. Therefore |/„(z)| > M for all z G C„. 

By taking N large enough we can ensure that for all n > A'^ that C„ is completely contained 
in U and that |^„(z)| < M/2 for all zinU. Therefore, for n> N and z G C„ we have 

\qn{z)\<M/2<\Uz)\ 

Hence by Rouche's theorem fn{z) and /,j(z) +^fi(z) have the same number of roots inside C„. 
Therefore p„(z) has a root Zn inside C„. Furthermore, Zn Zo because 

^ \Zn I ^ 

01 ,„ < argz„ < 02,« 

and rn,Rn 1 and 0i,„, 02,« ^ </>o- 

The remaining cases for when zo lies on one of the circles are very similar, but may involve 
an extra step. For example, if |zo — a| = 1 and |zo — b\ > I we apply the above analysis instead 
to the equation pn{z) /(z — = with 



/„ (z) = (z - a)" + j8 and ^„ (z) = a 



z — a 



+ 7 
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The final case is when zo lies on both circles, that is |zo — a\ = I and |zo — = 1. We 
first choose a sequence {Zq"'^ }m such that Zq'"'' ^ Zo with each Zo"*' belonging to only one of the 
circles. By the previous case, for each m there exists a sequence a sequence {zl'"' }« such that 
PnizT^) = and zi™^ ^ Zo""'. There exist positive integers Mi < M2 < M3 < . . . so that for each 
m we have 



I (m) (m) I 
\Zn Zq I 



< — for « > Mm. 

(hi) 



The sequence {zn} defined by putting Zn=Zn when M^ < n < M^+i satisfies Pn{Zn) = and 
Zn ^ Zo, as required. □ 

If one of the coefficients a or j8 vanishes then the form of the spectrum changes. This is 
illustrated in Figures [U and |2] and proved in Lemma[3] 

Lemma 3. Assume that <x ^Q, 15=0, 77^0 and b = —a, and that 

Pn{z) = {z-aY{z + a)" + a{z-aY + y. 
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IfO <a < \ then the roots converge asn^o°to the union of two arcs within the circles |z — fl;| = 1 
and \z + a\ = 1 together with an arc in the fourth order curve \ {z — a){z + a)\ = 1. The ends of 
all three arcs are at ib/-v/ 1 — a^. However, if a > I then the roots of pn converge to the entire 
circle \z + a\ = 1 together with the closed fourth order curve \{z — a){z + a)\ = 1 contained in 
{z ■ Re(z) > 0}. 



3r 
2 



1 - 



- 
-1 - 



-2 
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Figure 1: The constant a in Pn{z) of Lemma [3] satisfies < a < 1. 

Proof. First suppose that < a < 1. Let a^ = I, a2 = cc, aj, = y and 

/i(z) = {z-a){z + a) 
f2{z) = {z-a) 
Mz) = 1 

so that 

Pniz) = aifi {zT + a2f2{zT+a,fi{z)\ 
We will prove that as n ^ oo the roots of p„ converge to the union of the following subsets of C, 
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which we will denote by L. 



Li,2 = {zGC:|/i(z)| = |/2(z)|>|/3(z)|} = {z:|z + a| = l, Re(z)<0} 

Li,3 = {zeC:|/i(z)| = |/3(z)|>|/2(z)|} = {z:|(z-a)(z + a)| = l,Re(z)>0} 

L2,3 = {zeC:|/2(z)| = |/3(z)|>|/i(z)|} = {z:|z-a| = l, Re(z)<0} 

Lead = {±?Vl-a2} 

where Lend is the set of endpoints of the arcs Li 2, ^1,3 and L2,3. Let Lg be the e-neighbourhood 
of L. We will show that C\Le is contained in the union of finitely many sets of the form 

5,5 = {z G C : |/,(z)(l - 5)1 > max{|/,(z)| -.s^r}. 

Since Lg is bounded and since f\ (z) dominates both /2(z) and (z) in absolute value as |z| °°, 
there exists /? > and 1 > 5i > such that Lg c B{0,R) and z ^ B(0,R) implies 

|/i(z)(l-5i)|>max{|/2(z)|,|/3(z)|}. 

Therefore C\B{0,R) C BiSr Let K = B{0,R)\Le. For each zeK there is an G {1,2,3} and 
1 > 5^ > such that z G B^^ g^. Since each B^^ g^ is open and K is compact, we can cover K with 
finitely many such B^^ g^. In fact, if we take 5 small enough we conclude that 

3 

C\LgC |JB,5. 

r=l 

Inside the set B^g the term arfr{z)'^ dominates the remaining terms of Pn{z) so that for large 
enough n we have Pn{z) 7^ for all z G More precisely, if z G B^g then fr{z) / and 




> l-£K/a,|(l-5r 

> 

for n large enough, where we note that the bound is uniform for z G Bj.g. Consequently /?„ (z) / 
for all z G S^g for n large enough. Since C\Le is contained in the union of finitely many such 
Bj. g all the zeros of p„(z) are in Lg for n large enough, as required. 

We now show that if zo is in L then there is a sequence {zn} such that Pn{zn) = and Zn — zo- 
There are number of cases to consider, but the proofs are similar. We prove the case zo G Li^3 in 
a manner which emphasises the general technique. Since zo G Li^3 we have 

|/l(Z0)| = |/3(Z0)|>|/2(Z0)| 

By continuity there is an open neighbourhood U of zo such that |/2(z)//3(z)| < c < 1 for all 
z&U. We want to solve Pn{z) = for z close to zo, that is 

ai/i (z)" + 02/2 (z)" + 03/3 {zT = 0. 
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After dividing both sides by a\ and f3{z)" this is equivalent to solving 

(Mz)y a, a2(f2{z)Y 
\Mz)J ^a,^a,\Mz)J 

Putting /(z) =fi{z)/f3{z), a = a3/ai and^„(z) = (a2/ai)(/2(z)//3(z))" we want to solve 

f{zr + a + gn{z)=0 (2) 

for z near zo- We note that f'{zo) 7^ 0, f{zo) = e'^° for some 60 and gn{z) is uniformly exponen- 
tially small in U because |g,j(z)| < \a2/a\\c" where c < 1. 

Since f'{zo) 7^ the inverse mapping theorem implies that there exist open sets V and W 
such that Zo £ V, /(zo) G W and f :V is conformal with analytic inverse g -.W ^V. By 
reducing V if necessary we may assume that V <^U. 

Under the change of variable w = f{z) equation Q becomes 

w'' + a + gniw) = 

where gn{w) = gn{g{w)). We want to solve this equation for w near /(zo)- Since gn{w) is 
uniformly exponentially small in W as « ^ 00 we expect the solutions to be close to those of 

w" + a = (3) 

by Rouche's theorem. Let w„ = \a\^^"e"l'" be a solution of (|3]) with argument closest to 0o- 
This ensures w„ — > /(zo) as n ^ 00. Let C„ be the contour Cr„,i?„,e, „,e2„ with r„ = (^|fl:|)^/", 
/?„ = (||a|)'/", = 0,1 — 1^ and 02,;i = 0« + 5^- We note that w,, lies in the interior of C„. 
Routine estimates show that w G C„ implies 

\w" + a\ > \a\/2 

for all n. For « large enough we have C„ and its interior completely contained in W. Since 
gn (w) — > uniformly as « ^ 00 for w G VK, for « large enough we have 

\gn{w)\ < \a\/2 < \w" + a\ 

for all w G C„. Therefore by Rouche's theorem w" + a + gn{w) = has a solution m„ inside C„. 
Let Zn = gin,,)- The sequence {z,,} satisfies Pn{Zn) = and Zn ^ Zo, as required. 

The case a > 1 is handled very similarly, with the only difference being the set L which is 
now the union of the arcs 

Li,2 = {zGC:|/i(z)| = |/2(z)|>|/3(z)|} = {z:|z + a| = l} 

Li,3 = {zGC:|/i(z)| = |/3(z)|>|/2(z)|} = {z:|(z-«)(z + «)| = l,Re(z)>0} 

L2,3 = {zGC:|/2(z)| = |/3(z)|>|/i(z)|}=0 

^end = 

where we note that Li 3 is a closed fourth order curve which is actually completely contained in 
the open disk B{a, 1). The proof now continues as before. □ 
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Figure 2: The constant a in p„(z) of Lemma |3] satisfies a > 1. 

Lemma|4]is not covered by the general analysis and may be viewed as a degenerate case. 

Lemma 4. Let p be a polynomial in two variables of the form 

p{z,w) = amiz)w'" + am-iiz)w"''^^ -\ \-ai{z)w + ao{z). 

Let e > 0. As n ^ oo the solutions ofp{z,z") = that satisfy |z| < 1 — £ converge to the solutions 
of ao{z.) = satisfying |z| < 1 — e, while the solutions of p{z,z") = that satisfy |z| > 1 + £ 
converge to the solutions of am{z:) = satisfying |z| > 1 + £. 

Assuming that neither ao(z) = nor am{z) = has a solution of modulus 1, the remaining 
solutions of p{z,z") = are asymptotically uniformly distributed around the circle \z\ = 1 in the 
sense that for all a and j3 such that p — a <27l and for all 5 > we have 

#{z:p{z,z")=0,\-5<\z\<l + d,a<argz<l5} j8-a 

— ^ -r ' \ ; > as « ^ OO 

#{z:/7(z,z«) =0} 271 
where we count the roots with multiplicity. 

Proof We note that p{z,z") has nm + d roots counted with multiplicity, where d is the degree 
of am{z)- We begin by proving the last statement. Fix an arbitrary point zo of modulus 1. 
We first prove that the result is true in a neighborhood of z.q- The equation p{zo,w) =0 has 
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exactly m solutions counted with multiplicity because ^^(zo) ^ 0, and each solution is non zero 
because p{zo,0) = 00(^0) / 0- Let wi, . . . ,Wjt be the distinct solutions and let mi, . . . ,myt be the 

corresponding multiplicities so that m\-\ \-mk = m. The nm solutions of p{zo,z") = are 

precisely the rfi^ roots of the w,-, that is 

{z:pizo,f)=0} = \J{z:z"=Wi} 

i=i 

Let Zi be an w'*^ root of w,. We show that the solution Zi of p{zo,z") = has multiplicity m,-. 
Substituting z" for w and z" for w,- in 

p{zo,w) =a„,{w-Wi)'"'+a„+i{w-WiP+^ + ... {a„,^0) 

yields 

pizo,z'') = a„Xz"-zir+a„+i{z"-zir+' + ... 

= {z- ZiT' iz"~ ' + z"-hi + ■■■+ zr'r (am. + a„+i (z" - z?) + • • • ) 
= {z-Zirg{z) 

where g{z) is analytic and satisfies g(z,) / 0, and so the root zt has multiplicity m;. 

Let Ci , . . . ,Q be contours of the form Q^.e, .$2 such that the regions enclosed by the Q are 
disjoint and such that w, lies inside the region enclosed by Q. The preimage of C,- under the map 
z ^ z" consists of n disjoint contours of the form Qs^ei.ej which are equally spaced around the 
unit circle in the sense that they may be labelled c/'", . . . ,C"'" so that 

C/'" = C^-^C/'" for j = l,...,n where ^ = e'^''/" 

where multiplication by ^ effectively rotates the contour by an angle of iTi/n. As n ^ °o the 
radii of the preimage contours converge to 1 and their angular sweep converge to 0. We may 
label the roots of w; as z-'", . . . ,z"'" so that z.j'" is inside the contour C/'". 

Let C = U^^jQ. If w G C then p{zo,w) / since the only zeros of this polynomial are the 
Wi. Since C is compact there is a constant Mi > such that |;?(zO)>^)| > Mi for all w eC. Note 
that if z G Cj'" then z" is on C and so \pizo,z")\ > Mi. Let 

M2 = max{|w*| :wGC, 0<A:<m}>0 

The continuity of the polynomials a,(z) imply there exists an open neighborhood of zo^ say 
A, such that 

Ml 

\ai{z) - a,(zo) I < XT — -TTTT for < ? < m and z e A 

z(m + 1 jM2 

By reducing A if necessary, we may assume that 

A = {z : s < \z\ < S, < argz < (h.} 

Inside the contour C/'" the equation p{zo,z") = has precisely m, solutions, and these are 
all at the point z/'". If the contour Cp" Ues completely inside the region A and if z € Cj'" then 
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\p{zo,z") - p{z,z")\ < K{zo)-aUzW"\ + --- + \ai{zo)-ai{zW\ + \ao{zo)-ao{z)\ 

Ml 



< (m+1) 

< Ml 

< Ipiz^z") 



2(m + l)M2 



M2 



Therefore by Rouche's theorem the equation p{z,z") = has precisely m, solutions inside 
the contour C/'". Since the angle between successive contours C/'" is Ik/u we conclude 

#{y:CrcA} = !*|^ + 0(l) 

as n ^ °o where the error term 0(1) is actually < 2 for all large enough n. Consequently the 
number of solutions of p{z,^) =0 that lie in A is at least 

E.,x#0:CrcA} = !!^^^+O(l) 

Since the unit circle {z : |z| = 1} is compact we may cover it with finitely many sets of the 
form A, say A\,... ,An, such that in each A,- we have the above lower bound on the number of 
solutions of p{z,z") = in A,-. After trimming these sets and relabeUng we may assume that for 
i = l,...,Af that 

A; = {z : 1 - 5o < |z| < 1 + 5o, 0,-1 < argz < ^i} 

where So < S and (^0 < 0i < • • • < ^'a' with 0a? — 0o = 2;r. For n large enough we conclude that 
the number of solutions of p{z,f) = in the aimulus {z : 1 — 5o < |z| < 1 + do} is at least 

L 27r + 0{l)=nm + 0{\) 

Since p{z,z") = has only nm + d solutions we conclude that all but 0(1) of these solutions 
Ues inside a contour C/'". Therefore, the number of solutions of p{z,^) = in a region T 
can be counted with increasing accuracy as n — >^ 0° by identifying those contours C/'" which 
are contained in T and recalUng that the region enclosed by C/'" contains exactly m,- solutions. 
Therefore 



#{z : ;7(z,z") = 0, 1 - g < |z| < 1 + g, « < argz < /3} 
#{z:p(z,z") = 0} 



i:|Linm,-(j3-«)/27r + 0(l) 

nm + d 
nm{P-a)/27t + 0{l) 



nm + d 



15-a 
2n 
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as « ^ oo, as required. 

We now return to the first statement. Let B = {z : \z\ < I — £} and bn{z) = p{z,z") — ao{z) so 
tliat p{z,z") = bn{z) +ao(z). Let zi, . . . ,2^; be tlie roots of fl'o(z) in B. Witliout loss of generality 
we may suppose that e is small enough so that all the roots of ao (z) with modulus less than 1 are 
in B. For / = 1 let D, be a disk of radius 5 and centre z, , and let 5, be its circular boundary. 
For 8 small enough each D,- C B. There exists M > such that |<3o(z)| > M for all z in B\ U, D;. 
Since bn (z) uniformly on 5 as « ^ oo, we have for large enough n that 

\ao{z)\ >M/2>\bn{z)\ 

for all z in B\ U/D,-. Consequently the only solutions of p{z,z") = in B are in U,,D,-. Further- 
more, Rouche's theorem applied to the contour St implies that each D, contains a solution of 
p{z-,z") = 0. Since 5 is arbitrary the result follows. The proof of the second statement is similar. 

□ 

Example 5. Consider the polynomial p given by 

p{z,w) = + ai{z)w + ao{z) 

where ai (z) = — z^ — z + 9/2 and ao(z) = z^ — z^/2 — 4z + 2. Letting z = e'^ and solving for w 
gives two distinct solutions 

for all B e [0,271]. Since 1/2 < \f+ie'^)\ < 3/2 and 3 < \f-{e'^)\ < 5 we should anticipate that 
for large n the roots of p{z,z"') = will form two distinct rings of solutions, both close to the 
unit circle, along with an extra zero close to 1/2. This is illustrated in Figure [3] 

Example 6. Consider the polynomial p given by 

p(z,w) = w^ — 4w — 8z + 3 
Letting z = e'^ and solving for w gives 

w = 2 ± \/8e'« + 1 

For each 6 let w\ {6) and W2{6) denote the two roots. The function / : [0,27r] — > R defined by 

/(e) = ||wi(0)|-|w2(0)|| 

measures the difference in absolute value of the two roots when z = e'^ , and is well defined 
despite their being no canonical choice for wi (0) and W2(0). For each 6 there are two (possibly 
overlapping) rings of solutions of p{z,z") = in the vicinity of e'^ corresponding to those roots 
close to the solutions of z" = w\{6) and z" = W2{6)- If /(^i) > /(^2) the two rings of solutions 
near e'^' are further apart then when near e'^^ for fixed n. A straightforward calculation shows 
that 

/(e) =2(16cos + 65)1/^ 
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Figure 3: Solutions of p(z,z") = for n = 40 in Example |5l 



and so / decreases from to tt and increases from n to 271. Therefore one would expect the two 
rings of solutions to be closest together at — 1 , with the distance between the two rings gradually 
increasing as one moves from —1 to 1 along the rings. This is illustrated in Figure |4]for n = 20. 

The correct context for understanding Lemmas |2] and [3] involves functions of the type 

m 

Fn{z) = Y^ar{z)fr{zT (4) 

r=\ 

where {/r}"Li ^rid {arYr=\ '^^'^ analytic functions on an open connected set R. We 
impose the following hypotheses on the f/. 

1. There exists a compact set A" C /? and constant c € (0, 1) such that z G R\K imphes 

|/.(z)//i(z)|<c ifr^l. 

2. For each r ^ s the function fr/ fs is not constant. 

3. For all choices of the three distinct integers r,s,t the set 

{zeR:\fr{z)\ = \m\ = \m\} 

is finite. 
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Figure 4: Solutions of p{z,z") = for n = 20 in Example |6l 



Hypothesis 3 is false for the adjacency matrices of ordered graphs. In this degenerate case 
the functions fr are all powers of z, and there is only one curve involved, the circle |z| = 1, as 
seen in Lemma HI 

The set L defined in our next theorem is the union of what might be called the anti-Stokes 
lines of the function F„. We note that Hypothesis 3 is not required for the proof and so the result 
holds in a more a general context. 

Theorem 7. Let L (^R be the compact set consisting of all points at which two of the functions 
f. have the same absolute value and no others have a larger absolute value, that is 

m 

L=[j{zeR: \fr{z)\ = max{|/,(z)| : s + r}}. 

r=l 

Then given e > every zero of Fn is contained in the e-neighbourhood Lg of L provided n is 
large enough, with the exception of certain zeros which converge to certain zeros of one of the 
functions ar as n ^ oo. 

Proof. The set L is bounded because \ f\ {z)\ > \ fr{z)\ for r 7^ 1 and z outside a compact set, by 
hypothesis. It is also closed because each function \friz) \ — max{|/i(z)| : * 7^ r} is continuous in 
z. Hence L is compact. The set R\Le is contained in a finite union of sets of the form 

Br,S = {zeR: |/,(z)(l - 5)1 > max{|/,(z)| : s ^ r}} 
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where 1 > 5 > 0. In fact, by taking 5 small enough we conclude that 



r=l 



In the set B^g the term ar{z)fr{z)" dominates all the others for large enough n so / 0, 
unless z is close to a zero of a^. We make this precise as follows. Let 

m 

Z=\J{zeR:ariz)=0} 

r=\ 

and let Zg/ be the e'-neighbourhood of Z, so that in fact Zg/ is the union of the open disks B(zo,£') 
where zo is a root of an a^(z). Let A' be a compact subset of R and let C be the compact set K\Zgi. 
Let 



M = max{sup 

zee 



ai{z) 



and choose A'^ large enough so that 



{m-l)M{l-5f <-. 



If z G C\Le then there exists r G {1 , . . . ,m} such that z G B^g. This implies 



ar{z)fr{z} 



for n> N. Therefore for n large enough the only solutions of F„(z) = in K\Lg^ are within a 
distance s' from a root of one of the a^. Enlarging A' and decreasing e' shows that the zeros of 
F„{z) not in Lg converge to the zeros of one of the functions Ur as n ^ °o. □ 

We note that the zeros of that do not lie in B^ g for some 5 > have no spectral significance 
as « ^ oo. If however ar{zo) = and zo G B^g for some 5 > then there is a seqeunce {z„} such 
that Fn(z„) = and z„ zo- 

A function / : (m, v) — >^ C is said to be real analytic if it is C°° and if for all a in {u, v) there is 
a corresponding 5^ > such that 

is an absolutely convergent series for |;ic| < da. Note that / is real analytic if and only if Re / and 
Im/ are real analytic, and that / is real analytic if and only if / is the restriction of a complex 
analytic function f:U^Cto (m,v), where U is some complex neighbourhood of (m, v). Sums, 
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products, quotients and compositions of real analytic functions are real analytic, and the square 
root of a strictly positive real analytic function is real analytic because log has an analytic branch 
defined on C\[— oo,0]. 

The next theorem illuminates the structure of L. The following subsets of R will play an 
important role. We put A = Ai U A2 U A3 where 

Ai = {z G /? : {fr/fs)'{z) = for some distinct r, s} 

^2 = {z G /? : ar(z) = for some integer r} 

A3 = {zeR: \fr{z) I = \fs{z) I = \f, (z) I for some distinct r, s, t} 

Theorems. The set 

m 

L=[j{zeR: \fr{z)\ = max{|/,(z)| -.s^r}} 

is the disjoint union of sets D and Cj, j = \,. .. ,N, with the following properties. D is a finite 
subset ofC. Each Cj is a real analytic curve, i.e. the range of a real-analytic one-one function 
Yj : {uj,Vj) —> C whose derivative does not vanish anywhere. These curves satisfy the following 
properties: 

1. Yj is contained in 

{zeR:\fr{z)\ = \fs{z)\>\f{z)\} 
for distinct integers r, s (dependent on j) and all t ^ r, s. Ifz G Yj then {fr/fsYiz) 7^ 0. 

2. Yj ^'^y be completed to a closed curve, or its ends lie in D n A. 

3. Yjn{A2UA3)=idforj=l,...,N. 

Proof. We begin by showing that the result is true for each set 

L, = {zeR: \fr{z)\ = max{|/,(z)| : s / r}}. 

The first step is to show that for each zo in there is a corresponding open neighbourhood 
such that Lr n f/zo\{^o} is contained in the range of finitely many analytic arcs each of which has 
range contained in Ly. There are a number of cases to consider. 
If Zo ^ ^3 then there exists a unique s not equal to r such that 

\fr{Z0)\ = \fs{Z0)\ > \ft{Z0)\ 

for all t / r, .v. By continuity there exists an open disk centered at zo, say B^, and constant c > 0, 
such that z G B^^ implies 

\fr{z)\>C>\f{z)\ 

for all t ^r,s. Consequently the set LrDB^ consists of precisely those points z in B^ such that 
= or equivalently |/.(z)//.(z)| = 1. Define / : ^ C by /(z) = fr{z)/fs{z). 
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Since |/(zo)| = 1 there exists 0o £ R such that f{zo) = e'^- There are now two cases to consider 
depending on whether f'{zo) is zero or non-zero. 

If f'{zo) 7^ then by the inverse mapping theorem there exist open neighbourhoods U of zo 
and V of f{zo) so that f -.U has an analytic inverse g:V ^U. We may assume that y is a 
small disk of radius r < 1/2 centered at f{zo) and that U C B^^. There exist m, v G R satisfying 
u <v <u + 27t such that 

vns^ = {e^ ■.te{u,v)} 

where denotes the unit circle centered at 0. Define 7: (m,v) ^ C by /(f) = g{e"). We note 
that 7 is a real analytic arc because it is the restriction of g{e'^) to (m, v). Since 

LrCiU = {7(f) : t G (m,v)} 

it is sufficient to define U^g to be J7. 

If f'{zo) = things are a bit trickier. Suppose 

/(z) = ao + a^(z - zoT + «m+i {z - zq^^^ + ... 

where m>2 and 7^ 0. Obviously ao = f{zo) = e'^- We show that f{z) essentially behaves 
as z<-^ for z close to zo- We now make this precise. Let 

k{z) =am + flm+l (z - zo) + am+2(z - Zof + ■■■ 

Since k{zo) / we see that in a small enough open neighbourhood of zq, say U^, we can take 
an analytic branch of k{zY/'^. Let r(z) be such a branch. Now 

f{z) = aQ + {z-zo)'"{am + am+i{z-zo)+am+2iz-zo)^ + -..) 
= ao + iz-zork{z) 
= ao+[{z-zo)r{z)r 
= ao + h{zr 

where h{z) = (z - Zo)''(z) satisfies /i(zo) = and /i'(zo) = ?'(zo) / 0. Since /i'(zo) / we may 
shrink i/j.^ so that U^^ C B^.^^ and h : J/^, ^ Uq is conformal, where G J/o = h{U^^^). Choose 
e < 1/2 so that S(0,e) C Uq, and shrink [/^p again so that f/^g = /j^'(S(0,£)). This enables us 
to factorise /(z) as the composition of the following three surjective maps 

i^S(0,e) i^s(o,e-) i^£±^fi(ao,e-) 

Let s : B(0,e'")\[0,— ao] — 5(0, e) be an analytic branch of z 1— *^ z^/*" where the branch cut 
[0, — ao] is the Une connecting to — ao. The m analytic branches of z 1— z^^*" on S(0, e'")\ [0, — ao] 
are then given by 

Sk[z)=e'" 5(z) for^= l,...,m. 
RecalUng that ao = e'^^, there exist 0_ , 0+ G R such that d-<dQ<d+<d-+ In and 

S(ao,e'")n5i = {e'' G (0-,0+)} 
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For k=l,...,m define ak : 0o) — C by 

Ukit) = h-\sk{e^ - ao)) 
Similarly, for k = I,. . . ,m define jS^t : (^0) 0+) ^ C by 

l5k{t)=h-\sk{e''-ao)) 
We note that a^t and Pk are real analytic curves that satisfy 

m 

Lr n U^\{zo} = U {a*(0 : t e {9., Bo)} U {Pkit) : ? G 0+)} 

and so n f/j-g minus the point zo is contained in the range of finitely many analytic arcs, as 
required. 

Now we consider the cases when zq € A3. Since zo € Lr, fr and at least two other functions, 
say for j = 1, . . . attain the maximum absolute value, i.e. 

\fr{zo)\ = \U (zo)| = • • • = \U{zo)\ = max{|/,(zo)| : ^ / r} 

There are two cases to consider depending on whether or not fr{z()) is 0. We first suppose 
fr{zo) / 0. By continuity there exists £,c > so that if z G B{zo,£) then 

|/.(z)|>c>|/,(z)| 

\fsAz)\>c>\f,{z)\ iori=\,...,k 

for t 7^ r,Ji,...,j,fc. If necessary reduce e so that B{zQ,£)r\A-i = {zq}. For each i = 1,...,A: 
there exists an open set Ui C B{zQ-,e) such that the set {z G ?7/\{zo} : \fr{z)\ = 1/^,(^)1} can be 
parameterised by finitely many analytic arcs, say y['\..., yli} . Let 7 : (m, v) ^ C be one of these 
curves so that \ fr{Y{^)) \ = l/s,(7(-^))l fo^" all;c G (m,v). We show that either /CL^ or /nL^ = 0. 
Pick? ^ r,Si and suppose that there exists xi,X2 G (m,v) such that 

m{xi))\<\fMxi))\ 

and 

|/.(y(x2))|>|/.(r(x2))|. 

By the intermediate value theorem there is a point X3 between xi and X2 such that 

|/.(7(^3))| = |/.(7(X3))| 

and so 7(^3 ) G A3. But by construction 7nA3 =0 and so this is impossible. Hence we either 
have 

\Mz)\<\Mz)\ forallzGy 

or we have 

> l/r(^)l for all zGy. 
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If there isat ^ r,Si such that the latter case holds then ynL^ = 0, otherwise 7 C Lr. Let 

k 

i=i 

and 

k 

It is clear that 

L,nf/^„\{zo}c U rcL, 

as required. 

The remaining case is when /r(zo) = 0- This implies all of the functions fs vanish at zo 

and so one of these functions has a zero of lowest order, say k. For s = I,... ,m put gs (z) = 
/v(z)/ (z - Zo)''. Not all the g,{z) vanish at zo- If \gr{z[))\ = max{|^,(zo)| : * / 4 then gr{zo) ^ 
and so we apply one of the previous cases to the gs since for all z we have 

\fr{z)\ = max{|/j(z)| -.Sy^r} if and only if \gr{z)\ = max{|gj(z)| : J 7^ r}. 

If equaUty does not hold, then there is an open neighbourhood of zo such that Lr n = {zo}- 
We have now shown that for each z in Lr there is an open set containing z such that 
Lr n J/^VIz} is contained in the union of finitely many analytic arcs. If z is not on one of these 
analytic arcs then z G A. Since Lr is compact it can be covered by finitely many such open sets, 

say i/^j , . . . , i/^^, and consequently Lr is the union of finitely many analytic arcs, say 71 , , 7„ 

and a finite subset of A. We note that each 7 has non-vanishing gradient since each is of the 
form {fr/ fs)^H.s") for some analytic branch of (fr/fs)^^- However, as it stands, the 7 might 
not be disjoint from themselves and this must now be rectified. 

Now, suppose that the 7 are not disjoint. Without loss of generality we may suppose that 
7i : (mi,vi) ^ C and 72 : (m2,V2) ^ C are not disjoint, that is, there exists ti G (mi,vi) and 
h G ("2,^2) such that Yi(ti) = Jiih) = wofor some wq G C. We know that Y\{t) = (/^/Zs, (e'') 
for some and that 72(;) = {fr/fs2)^^{^") for some S2. If ^'i / S2 then wq G A3 which contradicts 
7 n A3 = for all i. Therefore si=S2 = s. Now 

e'" = {fr/fs){wo)=e''^ 

and so ti = ?2 + 2;rn for some integer n. Since e" is 2k periodic we can reparameterise the range 
of 72 by 

72 : (M) ^ C where 72(0 = {fr/fsY'ie'') 

where 1*2 = U2 + 2nn, V2 = V2 + 2nn and where [fr/ fsY^ is the same analytic branch used in the 
definition of y2. This ensures that 7i(?i) = = wq. Let f = fr/fs and note that /(wq) = e''^ 
Since /'(wq) 7^ there exist open neighbourhoods U of wq and V of e"^ such that f :U has 
analytic inverse g.V—>U.IfS>Ois small enough then for all t in {ti — S ,ti + 5) v/e have 
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Since 71 and 72 are analytic and agree on a small interval we deduce that they must also agree 
on all of (mi, vi) n {u2,V2) by the principle of isolated zeros. Hence we can replace 71 and ^ by 
the single analytic arc 



T : (min{Mi,M2},max{vi,V2}) — C where t(?) 



7i(0 if/G(Mi,vi) 
72(f) otherwise 



If we continue the above process of replacing two overlapping arcs by a single analytic arc 
then after finitely many repetitions we will arrive at finitely many disjoint analytic arcs, say 
such that U,T, = UiYi. It is clear that for each t,- there is a corresponding j,- such that 
the range of T,- is contained in 

{zG/?:|/.(z)| = |/.,(z)|>|/,(z)|} 

for all t ^ r,Si. 

There is no guarantee that each T,- : (m, ,v,) ^ C is one-one, but this situation can be easily 
rectified. Suppose T; : (m,-, v,) ^ C is not one-one. There then exist a,be {ui,Vi) with a<b such 
that T,(a) = Ti{b) and Ti\{a,b) is one-one. Since 

e'" = ifr/fsm^)) = {fr/fs.){<b)) = 

we conclude that b — a = Inn for some n. Let f, : (m,,v,) ^ C be the unique 2nn periodic 
function such that Ti\[a,b) = Ti\[a,b). If 5 > is small enough then for |/| < 5 we have 

z{a + t) = 'c{b + t). 

This impUes T; is analytic and so T,(f) = f,(f) for all t G {ui,Vi). Therefore we may replace T, 
by the one-one function Zi\{a,b) and its endpoint T;(a). So we can assume that all the T, are 
one-one. 

We conclude that for each r = 1, . . . ,m the set is the disjoint union of a finite subset Dr 

(r) (r) 

of C and one-one analytic curves Yi ,---,Ynr that satisfy Properties 1 and 2 in the statement of 

(r) (s) 

the theorem. It is possible that for r ^ s that Yi is not disjoint from Yj for some i and j. It 

(s) (s) 

follows by Property 1 that 7- C L^. Hence y- is redundant and can be deleted. We continue 
this process until the remaining curves are disjoint. We conclude that L = U^L^ is the disjoint 
union of a finite subset Z) of C and one-one analytic curves, say 71 , . . . , 7„, that satisfy Properties 
1 and 2. 

Finally, we want each 7 to also satisfy Property 3, that is 7 n (A2 U A3) = 0. We first note 
Ln (A2UA3) is finite because L is compact and A2UA3 discrete. If 7 : (m,, v;) C meets A2U A3, 
that is 7 n (A2 U A3) / 0, then there exist finitely many points, say ai , . . . , a^, such that 

Ui<ai <a2<-- - <am< v, 

and Yi{t) € (^2 U A3) if and only ifte {ai ,am}. We now just replace Yi by 

ri\{ui,ai), Yi\{ai,a2),..., Y\{am,Vi) 

and enlarge D with the points Yii^k) for all k. After doing this for each Yi, we obtain a finite subset 
D of C and a finite collection of analytic curves that satisfy all the required properties. □ 
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We now show that each arc 7 : (m, v) ^ C in Lemma [8] can be parameterised by arc length 
with the arc length parameterization also being real analytic. We begin by showing that the 
length of 7 is finite, that is 

/= f\y{t)\dt <oo. 

This is achieved by showing that 1/(01 = 0{{t — u)") as t ^ u for some a > — 1, and likewise 
at the opposite endpoint v. The problematic case is when for t near u we have 

Y{t)=h-\sie''-e'")) 
where h is conformal and 5 is a branch of z 1-^ z'/"'. Differentiating it suffices to show that 

\s'{e'' -e'")\=0{it-uy') 

for some a > — 1. Now 

m 

< co(? 

for some constant cq > as f ^ m, as required. Hence Z < 00, as required. This implies the length 
function s : (m,v) — > (0,/) where 

^(0= f\yi^)\dx 

Ju 

is well defined. Since 

and Y{x) does not vanish we conclude that is real analytic. Therefore s is real analytic. 

Since s'{x) > for all x S (m,v) we conclude that s^^ is real analytic by the inverse mapping 
theorem. Therefore the arc length parameterization 7 : (0, Z) ^ C of 7defined by f{t) = y{s^^ (?)) 
is real analytic. 

Our next theorem analyses the distribution of the zeros of F„ along L. The proof identifies 
the exceptional points and the densities in the statement. 

Theorem 9. There exists a finite subset D of L, real-analytic curves : (0,/^) C for r = 
. ,N parameterized by arc length with non-vanishing gradients, and bounded real-analytic 
densities Pr : (0,/^) (0,°°) with the following properties: 

1. L is the disjoint union ofD and the ranges of all the 7^. 

2. For every choice ofr and every closed subinterval [ot, j8] of (0, Ir) we have 

hm Ne= Pr{s) ds, 
e^O Ja 

where 

Ns=lim -\{z:F„{z)=0}nSe\ 

n^oo fi 
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and 



5e = {zGC:dist(z,7,([a,j8])) <£}. 



Proof. Let 7 : (0, /) — > C be one of the analytic curves in L. It has range contained in 

{zeR:\Mz)\ = \Uz)\>\Mz)\} 

for some integers r s and all t ^r,s. By taking a sufficiently small open neighbourhood U of 
7([a, j8]) we can ensure that z^U implies ar{z) / 0, /s(z) / and 



< 1-5 



for some 5 G (0, 1) and all r,s. Since ar{z)fs{z)" / for all z G ?7 the solutions (counted 
with multiplicity) of Fn{z) = in ?7 are precisely the same as those of F„{z) / {ar{z)fs{z)") = 0. 
Now 

Fniz) _(Mz)Y asiz) y a,{z) (mV 
ar{z)Mzr \Mz)J ^a,{z) ^,fca,(z) \Mz) J 

and so with/(z) = fr{z) / fs{z), a{z) = as{z) / a^iz) ^"(^^ gn{z) = Zk^r,Mk{z) / ar{z)){fk{z) / fs{z))" 
we are interested in solving the equation 

/(zr+a(z)+g„(z)=o 

inside U . We note that g„(z) is uniformly exponentially small inU as « — > 00. 

Define the argument function : (0, /) — > R by d{t) = — /log(/(7(f)) using an analytic 
branch of log/o 7. We note that Q is real analytic and can be extended continuously to the 
closed interval [0,/]. Also, define the density function p : (0,/) ^ R by p{t) = 6'{t)/27t. Sup- 
pose 7is the arc length parameterization of T : (m,v) C where T is of the form f^\e") as in 
TheoremlU Letting s : (m,v) —>■ (0,Z) be the length function 

s(t) = ( W(x) \ dx 



we recall that 7(f) = T{s^^{t)). Consequenty d{t) = s^^{t) +c for some constant c and so p{t) 
is bounded if {s^^)'{t) is bounded. It suffices to check that there exists a constant C > such 
that |/(?)| > C for all t £ (m,v). Since 

\s'{t)\ = \z'{t)\ = \{f-'ne')\ 

where is an analytic branch of a local inverse of /, it suffices to check that /'(z) is bounded 
for z G T (or equivalently 7.) But if one observes the construction of T this is clearly the case. 
Hence p is bounded. Also p is positive because {s^^)'{t) is non zero and s is increasing. 

The idea of the proof is to now cover 7([a, j3]) with finitely many sets with disjoint interiors 
such that in each set we can count the number of solutions of F„{z) = 0. The boundaries of 
these sets contribute 0(1) zeros as « — > 00 and so their contribution can be ignored. This will be 
achieved over a number of steps. 



21 



Let zo S 7([oJ)iS]) and let wo = /(zo)- Since f'{zo) ^ there are open neighbourhoods V of 
Zo and W of wq such that f :V —^W is conformal with analytic inverse g :W —^V. We may 
assume that V is sufficiently small so that V and z G V implies 

|a(z)-a(zo)| < 

Since |wo| = 1 there exists 00 G R such that wq = e'^°. Let C„,(, be the contour Cr.RSuBi where 
r = 1 — 5, /? = 1 + 5, 01 = 00 — 5, 02 = 00 + 5 where 5 > is small enough so that C^^^ and its 
interior is contained in W . Let A„,q be the union of C^^^ and its interior, so that A„^^ is a sector of 
an annulus with angular sweep 62 — 61- Finally let N^^ = g{Awg). We show that 



#{zEA.°:F„(z)=0} = ^^(^ + O(l) 



as « ^ 00. Let Fn{w) = Fn{g(w)) for w G VK. Since ^ : ^ V is conformal Lemma [T] implies 

#{z G A^; : Fniz) = 0} = #{w G : = 0}. 

The solutions of f„(w) = counted with multiplicity are precisely the same as those of 

w" + a{w)+g~n{w)=0 

where a(w) = a{g{w)) and gn{w) = gn{g{w)). We use Rouche's theorem to show that each such 

w" + a(wo) =0 



solution in A° is close to a root of 



which has roots 

wj = \a{wo)\^^"e''i forj=\,...,n 

where sj+i — sj = 2n/n. For j = I,... ,n define C^" to be the contour Cr,,^R,,^ei j,62j where r„ = 
(i|a(wo)|)i/", Rn = (||a(wo)|)^/", 6ij = sj - n/2n and 02j = sj + n/2n. We note that the O^" 
approach the unit circle uniformly as « ^ 00 and are evenly distributed around the unit circle by 
an angle of In/n. 

We show that if C^ " C A„g then Fn{w) = has precisely one solution inside C^ ". Routine 
estimates show that if w G C-''" then |w" + a(H'o)| > |a(wo)|/2, where we note that the lower 
bound is independent of n. Also, since gn{w) uniformly as « ^ 00 there exists N > such 
that for n> N we have |^„(w)| < |a(wo)|/8 for all w G Therefore for w G C^'" and n> N we 
have 



w" + a{w) + gn{w)\ > \w" + a{wo)\-\d{wo)-d{w)\-\gn{w)\ 
^ \a{wQ)\ \aiwQ)\ \diwQ)\ 
- 1 8 8 

_ |a(wo)| 
4 
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and 

\{w" + a{w)+gn{w))-{w" + a{wo))\ < \d{w) - a{wQ)\ + \g„{w)\ 

|a(wo)| \d{wo)\ 
8 8 
^ \d{wo)\ 
4 

Therefore by Rouche's theorem w" + a(w) +gn{w) = and w" + a(wo) = have the same num- 
ber of solutions inside C^'", and so F„{w) = has precisely one solution counted with multiplicity 
inside C^' ". 

Let D-i'" be the open region enclosed by C^ ". Routine estimates show that if w G A„g\UjD->'" 
then \w" + d{wo)\ > |a(wo)|/2. Consequently for w G Awg\ UjOJ'" and n> N we have 

K + a(w) + g„(w)|>^ 

and so F„{w) = has no solutions outside the C^'". Therefore 

#{zGiV° :F„(z)=0} = #{wGA°„:F„H=0} 
= #{y:C^>CA,J + 0(l) 

as « — > oo. We also note that the boundary of N^^ has 0(1) (in fact, at most 2) solutions as « ^ oo. 
Our choice of (and hence N^^) was quite arbitrary and it will be useful to impose further 
conditions. There exists to G [«, j3] such that zo = Y{to)- We may choose A^q so that there exist 
a, b satisfying 0<a<?o<^<^so that 

S'nA^, = {e''^<^'^ :te[a,b]} 

and 

rr\N,o = r{[a,b]). 

This uses the fact that f{y{t)) = e'^^'^ by definition. These sets satisfy the relation 

f{YnNj = s'nA,,. 

Furthermore we can impose the condition that d{b) — d(a) < n/2. 

For convenience, let (0i , 62) denote the closed region enclosed by the contour Ci_s^i+s,ei ,&2 ■ 
We show that there is a 5 > small enough and a partition of [a, J3] of the form 

to <ti = a <t2<t3 < ■■■ <tm = P < tm+1 

SO that /^^ has an analytic branch gi in a neighbourhood of A; = Ag{d{ti),d{ti+i)) so that the 
sets Ni = gi{Ai) for / = 0, . . . ,m have the following properties: 



23 



1. The interior of U,A^,- is an open set containing y([a,j3]). The interiors of each A^,- are 
disjoint. 

2. For j = 0, . . . ,m we have 5^ nA,- = {e'^W : t G [tuh+i]} and yn A^,- = 7([?i,?,+i]) 

3. For j = 0,...,m 

#{z E ^° : = 0} = + 

as n — oo. The boundary of A^, has 0(1) solutions as n — oo. 

Let z G 7([o:, j8]). By the previous step there exist subsets of C, denoted by A'^ and A^, and 
real numbers < b^, such that z G A^°, A^ = As^{0{az).,0{bz)) for some 5^ > 0, 

S' nA, = : t G [a„Z7,]}, yn^, = 7([«z,^z]) 

and 0(^7) — d{az) < n/2. There are also open subsets of C, say and V^, such that C J/,, 
A, C V, and f -.Uy^V^ has analytic inverse : ^ i/,. Since 7([o:,j3]) is compact there exist 
finitely many points, say zi , . . . ,zm, such that 

M 

7([a,/3])cUA^° 

Let P be the set formed by removing from 

:/ = l,...,M}u{a,j8} 

all elements less than a, except for the largest such one, and similarly, all elements greater than 
j3, except for the smallest such one. The set P can be labeled in increasing order so that 

tQ <h = a <t2<h < ■■■ <tm = ^ < tm+l- 

Initially let 5 be the minimum of the 5^.. For i = 0,. . . ,m there exists Zk (depending on /) such 
that [ti,ti+i] C [a^^,h,J. Letting g,- = we put A,;^ = As{e{ti),d{ti+i)) and Ni^s = gMi.s)- 

It now suffices to show that by reducing d if necessary we can ensure that the interiors of the 
A^,- 5 are disjoint. Suppose there exists / < j such that A^°g <^Nj g 7^ 0. This implies < 0{tj) 

for if 9{ti+i) = 9{tj) then we would have A°g nA^ ^ = (because 9{tj+i) — 9{ti) < n) and so 
N°gnNjg = 0, a contradiction. Thus and [tj,tj+i] are disjoint and hence 7([f,-,?i+i]) 

and Y{[tj,tj+i]) are disjoint. Consequently there exist disjoint open sets J7, and Uj such that 
7{[ti,ti+i]) C Ui and Y{[tj,tj+i]) C Uj. Now reduce d so that gi{Ai^s) ^ Ui and gjiAj^g) C Uj. 
This implies A^; ^ ^^j,S = ®- We continue this process of reducing 5 until the interiors of the 
A^,- 5 are disjoint. We then put A, = A, g = Ag (0 (?,) , 9 (?;+i )) and A^; = g,(A,) and it is clear that 
these sets satisfy the required properties. 

By definition Sg is the £ -neighbourhood of 7([o;,j3]). There exists c > such that for each 
j = 1 , . . . ,m — 1 we have 

f{NinSe)^Ac{9{ti),9{ti+i)). 
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Therefore as n — ^ 0° 

Hz^N?nS.:F^{z)=0} = <^^^±p^ + O{l). 

Also, for each £ > small enough there is a corresponding Cg > such that Cg — > as e — > 
and such that 

f{NonSe)CAc{e{h)-Ce,e{h)) 

and 

f{N„,nSe))CAa{e{t„),e{t„) + Ce) 
Therefore as n — >^ 0° we have for j = 0,m 

Ki,e := #{z eN?nSe:F„{z)=0}<'^ + O{l) 

Consequently, asn^°o 

m 

|{z:F„(z)=0}n5e| = £#{zGA^°nSe:F„(z)=0} + O(l) 



Hence 



and so 



Finally, 



and so 



e(^)-e(«) ^ ^,^^ e(/3)-e(«) ^ c, 

271 ~ ~ 271 71 

e^O 27t 

Pd'{s)^ 6(15) -d{a) 



27t 



lim Ne= p(s) ds 
e^O Ja 

as required. □ 

The above result is sharp in the sense that one cannot let a = or j3 = Z^. To see this consider 
the function 

Fn{z)=h{zr+f2{zr+Hzr 

where /i (z) =z2(^- 1), f2{z)=z- 1 and/3(z) = (z- 1)^. Since |/i(l)| = |/2(1)| = |/3(1)| =0 
we note that 1 does not lie on any of the analytic arcs. The problem is that = has n 
solutions counted with multipUcity at 1. If z is close to 1 and satisfies |z| = 1 we have 

l/l(^)l = 1/2(^)1 >I/3(Z)I 
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and so there is an analytic arc 7 : (0, /) ^ C with endpoint 1. Without loss of generality suppose 
that 7(0+) = 1. If we choose j3 > small enough we get 



This makes use of the fact that 6 can be extended continuously to [0,/]. For each e > the 
£-neighbourhood Se of 7((0,j3]) contains 1 and so Ne > 1. Consequently 



Corollary 10. The limit set of{z '. F„{z) =0} as 00 equals the union ofL and those zeros of 
a,.{z.) that belong to some B,.sfor some 5 > 0. 

Example 11. Given c > 0, the solutions of the equation (z^ — 1)" = c are given explicitly by 



where I <r <n. Putting f{z) = — 1, the solutions converge as « ^ 00 to L = {z : |/(z)| = 1}- 
As in Theorem [8] one may parameterize L by y{u) = +e'" where —71 < u < n. For this 
initial parameterization one has f{y{u)) = e'", 6{u) = u, p{u) = 1/ {2n) and j^j^p{u)du = 1. If 
instead one parametrizes the curve 7 by arc length using the formula 



then one sees that the density of the zeros vanishes at z = because s'{u) ^ 00 as m ^ ±71. 

3 Spectrum of a Directed Graph 

In this section we introduce some relevant notions from graph theory and use them to study the 
adjacency matrix of a directed graph. As well as providing a simple introduction to the general 
theory, this case has some features of its own that are of interest. 

Let ^ denote the class of finite, directed graphs (5, -^) that are irreducible in the sense that 
for every m, v G 5 there exists a path u = s\ ^ S2 ^ ■■■ ^ Sn = v. The irreducibiUty assumption 
implies that the indegree and outdegree of every vertex of S is at least 1 . 

The adjacency matrix A of (5, — >) G §^ is defined by 



The spectrum of (5, — >) is by definition the spectrum of A. 

We define C to be the set of all ^ G 5 that have total degree 2. If a G C and b ^ a ^ e then 
irreducibility implies that ft 7^ a 7^ e if #(S) > 1. If C = 5 then 5 is a cycle and 



liminf Ne > p{s) ds. 





,IH|-l/2 




Spec(A) = {f'2'^'''/":r=l,2,...,«} 
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where n = #(5). We henceforth assume that C and / = S\C are both non-empty. 

The set C can be written as the union of disjoint 'channels' Q, which we define as subsets T 
of S that can be identified with {b = 1,2, ...,e — l,e} in such a way that 

1 . every xeT has outdegree l;ifx<e then x — *■ x + 1 ; moreover e ^ e e J; 

2. every x G T has indegree 1 ; if x > 1 then x — 1 ^ x; moreover J 3 h ^ b. 

We write C = ULi Q where Q is a channel of length with ends Z?;, G 7. 

The structure of J may be quite complex, but every x G / has total degree at least 3. We may 
write J as the disjoint union of subsets that are coimected with respect to the undirected 

graph structure inherited from (5, —^). It follows that there are no directed edges joining different 
subsets Ji, which we call junctions. Irreducibility implies that one can pass from any junction to 
another, but only via intermediate channels. In some cases one may only be able to pass between 
different points of the same junction via external channels. 

We now lengthen each channel C,- by the factor n without changing its end-points bi, e,- or 
the junctions /, to produce a new graph {S'^"\^) G Sf, with a corresponding adjacency matrix 
aW. We will prove that for every e > the spectrum of 5^"^ is almost entirely confined to an 
aimulus of the form {z : 1 — e < |z| < 1 + £} if « is sufficiently large. There may however be 
a small number of eigenvalues far from the unit circle. In a later section we will show that the 
corresponding eigenvectors are concentrated around the set J. In general they are concentrated 
around a single junction, but this may not happen if the graph has certain symmetries. 

Our general analysis depends upon collapsing each chaimel Cr to a single point pr, so that 
we are left with the set 

S = JU{pr: 1 <r<h}. 

If /, J G 7 then we write / ^ j m S if and only if / j in S; we also include the extra edges 
br —>■ Pr and pr — > ir. A is defined to be the adjacency matrix of {S,—>) G Note that (5,^) 
and A depend neither on n nor on the lengths Cr of the original chaimels. 

Theorem 12. The eigenvalues of A^"^ coincide with the roots of the polynomial pencil D{z)—A 
on S, where D{z) is the diagonal matrix with entries 



Dijiz) 



Indeed 



for all z G C. 



Z ifi = j£J, 

z'^^' ifi = j = Prfor some r G {I, 
otherwise. 

det(z/-AW) = det(D(z) - A) 



.,h}, 



(5) 



Proof. It suffices to prove the following more general result which will also be of use later on. 
Let M be a matrix of the form 



C2 



M- 



\Bi B2 



Ch 
Bh 



El \ 
El 

Eh 
J 
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where Q is the Uk x rik matrix 



/ ak -1 



CCk 



\ 



CCk J 



with ak on the diagonal, —1 directly above the diagonal and elsewhere; J is an arbitrary mxm 
matrix; is an m x % matrix with everywhere except possibly for entry (4, 1) where it has 
the value j3^; and Ek is m rik xm matrix with everywhere except possibly for entry {n^Jk) 
where it has the value Sk- It is clear that the vertices of the graph (5^"^ , — can be labeled so that 
the corresponding matrix zl — A^"^ has the above form. 

On the set of matrices of the above form we may define a map M M hy defining the 
corresponding matrix M by 



M = 



rii, 



E'2 



where is the m x 1 column vector with value in position and elsewhere; E[ is the 
1 X m row vectors with value £k in position jk and elsewhere; and where the 1 x 1 matrix (a^*) 
replaces the matrix Q. It is clear that the vertices of the graph (5, — >) can be labeled so that the 
corresponding matrix D{z)—A has the above form. 

It now suffices to show that det(M) = det(M). We prove the result by induction on h which 
is the number of matrices Q along the diagonal of M. The result clearly holds for /j = and for 
any size matrix J. For the inductive case, let A'^ = ni H Vnn and let K be the matrix 



K = 



( C2 



V B2 



El \ 

Ch Eh 
Bh J J 



Expanding the determinant of M down the first column gives 



det(M) = aidet(Mi,i) + (-l 
= ai"'det(^) + (-l 



Pi det(M, 



1 



where det(Mi i) is evaluated by expanding successive minors down the first column and det{MN+ii,i) 
is evaluated by expanding successive minors across the first row. Similarly, 

det(M) = det(^) + (-l)''+'i+i^i(-l)''-i+^'+iei det(^,_i+,-,,,_i+,-J 
= <i det(^) + (-!)'■' +^''+1)31 ei det(J^;,_ 1+;, ,h-i+n ) 
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By the induction hypothesis 

det{K) = det{K) 

and so det(M) = det(M), as required. 

Corollary 13. There exists a polynomial p of the form 

p{z.,w) = ac{z)w'-' + ac-\{z)w'^^^^ H a\{z)w + ao{z) 

such that 



□ 



(6) 



Att{zI-AW)=p{z,f). 

Moreover c = #(C). All of the polynomials ar have degree at most #{J) and ac has degree equal 
to #(/). 

Proof. Theorem[T2limplies det(z— A'")) = det(D(z) —A), and after the change of variable w = z" 
we note that 



det(D(z)-A) 



w 



£2 



B\ B2 



E'2 



B'u zI-J 



Using the permutation definition of the determinant we see that the highest order power of w 
corresponds to the term in the determinant that includes all of the terms and the power of 

w is then ei H \-eh = c. The term in the determinant arising from the identity permutation 

shows that ac{z) has degree equal to #(/). It is clear that this is the maximum degree of the 
polynomials ar{z)- □ 

It is easy to construct numerical examples which illustrate the asymptotic spectral behaviour 
proved in Lemma HJ but we give a number of more general examples below. 



4 The General Model 

We now turn to the general model described in the introduction. Given a matrix A we construct 
the associated graph and define the channels and junctions as described before, except for the 
following variations. We allow every point x in a channel to have an extra edge which starts 
and ends at that point provided the diagonal matrix entry Av,jc is the same for all points in the 
channel. We also require that all the non-zero matrix entries Av,v associated with a particular 
channel, including the two that link it to external junctions, should be equal. Thus every channel 
C, is determined by three parameters, its length et, the diagonal entries a,, which may or may 
not vanish, and the non-zero, off-diagonal entries j8,- associated with edges that link successive 
points in the channel. 

Given a positive integer n we now construct a new matrix A^") which has the same graph 
structure and matrix entries as A, except that the length of each channel is ne, instead of e,. We 
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are interested in the asymptotic behaviour of the eigenvalues of A^") as n — > oo. The method of 
Theorem [T2l implies that 

p„{z) :=det(z/-AW) =^£a,(z)n((z-«,)/A)"'^' (7) 

s ies 

where )k = nil A"''' the sum is over all subsets s of {\,2, . . . ,h} and each a, is a polynomial of 
degree at most d = Ly=i ^{Jj)- If <7 = {1,2, . . . then is of degree d. The other Us may be 
of lower degree and may vanish identically. The zeros of Pn{z) coincide with those of 

f„(z)=£«.(z)/.(zr 

s 

where 

/s(z)=n((2-«o/A-r' 

and so the spectral asymptotics are determined by Theorem |7] 

Since there are 2'^ subsets of {1,2, . . . one would expect the number of terms in Pn{z) 
to grow exponentially as the number of channels is increased. Surprisingly this often does not 
occur as many of the a, may vanish identically. We will prove this for the case when each 
junction consists of a single point. We will need the following coefficient theorem for directed 
graphs; see [6, Theorem 1.2]. Because of the presence of the term z/ we allow / i to be an 
edge for all / G S whether or not it was previously. A single point of S is regarded as a cycle of 
length 1. 

Theorem 14. If s C {1^2, ... ,h} and a,,- is non-zero then after removing every channel Ci such 
that i € sfrom the graph it is possible to cover the remainder by disjoint cycles such that each 
channel C not in s is completely contained in one of the cycles. 

Proof. We use the formula 

It ieS 

where A^ = zl —A and the sum is over all permutations n. We note that the term tji = Ylies^z.i,7i{i) 
is non zero if and only if (/,7r(/)) is a directed edge of S for every /. Every n may be written 
as a product of disjoint cycles, some of which may be of length 1. If one point of a channel 
C is a cycle of length 1 then the same applies to every point of C and tj^ contains the factor 
((z — a,)/j8,)"'. Otherwise the contribution of the channel C to tj^ is a non-zero constant. The 
polynomial as is therefore obtained by summing over those partitions of S into disjoint cycles 
such that every point of every channel in 5 is a singleton cycle, and each channel not in s is a 
part of a cycle of length greater than 1. Since a., is non-zero such a partition actually exists. □ 

Corollary 15. Suppose each junction in S consists of a single point. Let d be the number of 
junctions and let h be the number of channels. IfN is the number of subsets s of {I, . . . ,h} such 
that as is non zero then N = 0{h^) ash^ °°. 
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Proof. Let G., denote the graph obtained from S by removing each channel Q for / G s. If Gs 
can be covered by disjoint cycles such that each channel C not in s is contained in a single cycle 
then for each junction in Gs we either have precisely one in channel and one out charmel, or else 
we have no chaimels at all. So we are now left with a counting problem. The number of subsets 
s such that every junction in Gs has precisely one in channel and one out channel is at most 
(^) since each possibility requires precisely d channels, and we have h from which to choose. 
Similarly, the number of subsets s such that A: junctions in Gs have precisely one in chaimel and 
one out chaimel with the remaining junctions having no channels is at most (^) . Consequently, 

as ^ — >^ oo, as required. □ 

Despite the above result, there are examples for which Us is non zero for all 2^" subsets s of 
{1 , . . . , /i}. Let 5 be a graph consisting of one junction and h chaimels that satisfy the following 
properties. The junction 7 is a cycle of length 2h with 

1^2^3^ ^2h-\^2h^\ 

For j = 1 , . . . , ^ the channel C, begins at i and ends at i + 1 so that 

i — ^ C, — ^ j + 1 

The length of each channel is not important. If j is a subset of { 1 , . . . , /i} then the only way of 
covering the graph Gj by disjoint cycles so that each chaimel C not in j is contained in a single 
cycle is by using the unique cycle in Gs that contains all the channels. Therefore each is 
actually a non zero constant. 



5 Some Examples 

The matrices A '^"^ considered in this paper may be classified according to the number of channels 
h and the number of junctions k. In this section we consider some of the phenomena that arise 
for small values of h and k, our goal being to describe the asymptotic form of Spec(A("') as 
n — >^ oo. Since each junction has at least one in channel and one out chaimel by irreducibiUty, it 
follows that k <h. Even if we assume for simplicity that each junction consists of a single point 
and that each channel has the same length n, there are still several graphs for each choice of h 
and k. 

Example 16 {h = k= 1). We suppose that the graph has one junction containing two points and 
one chaimel containing n — 2 points. The nxn matrix is taken to have the form 



A = 



' 1 


if J = j+1, 


a 


if / = ;•= 1, 


b 


if / = 1 and j = n 


c 


if y = 1 and / = n 


d 


if / = = n. 


. 


otherwise. 
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in which the constants a,b,c,d describe boundary conditions at the ends of the ordered interval 
{1,2, . ..,«}. One may check that the characteristic polynomial of A is 

Pn{z)=f~^ {z^-{a + d)z + ad-bc)-c. 

The roots of this polynomial converge asymptotically to the unit circle, with the exception of 
isolated roots that converge to any solution of — (a + d)z + ad — bc = that satisfies |z| > 1 . 

Example 17 (h = 2, k = I). We consider the case in which the graph associated to A has only 
one junction and that junction contains only one point. We also assume that there are h = 2 
channels, each of which has length n. The matrix A therefore has 2/2 + 1=5 free parameters 
apart from n. We will see that all the anti-Stokes curves are circles, where we regard straight 
fines as circles with infinite radii. For n = 3 the matrix A is of the form 

\ 





/ "1 


Pi 


Pi 

CCl 








Pi 


A = 










Pi 


Pi 


Pi 






Pi 


r 



where we need to assume that j3i and j32 are non-zero. For general n the reduced matrix is of the 
form 

/ (z-ai)" -j8f 
A{z)={ (z-a2r -PS 

\ Pi Pi z-r 

and the characteristic polynomial is 

Piz) = {z-r){z~ aiYiz - - ft"+' (z - «2)" - Pt' (z - ai)". 

This may be simplified further if ai = a2 and we assume that this is not the case. 
The characteristic polynomial is of the form T,r'^r{z)fr{z)" where 



fiiz) 
fiiz) 

Hz) 



(z-ai)(z- 

jSl(z-«2) 

Pi{z-a,) 



and a\{z) = z — y, ciiiz) = —Pi, a'3(z) = —Pi- Each of the functions fi dominates the others in 
absolute value in an open set Ui, where Ui contains all large enough z, cci € U2 and a2 € U3. 
The limit set E of the spectrum of A as « ^ 00 is contained in the union of the anti-Stokes lines, 
with the possible exception of an eigenvalue that converges to 7. Whether or not this eigenvalue 
exists depends on the parameters of the matrix. 
The anti-Stokes lines are defined for / ^ j by 



Kij = {z:\fi{z)\ = \fjm 
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and are given by the formulae 



^1,2 = {z:|z-ai| = |j8i|} 
^1,3 = {z:|z-a2| = |i82|} 
^2,3 = {z:|i8i||z-a2| = |j82||z-ai|} 

An elementary calculation shows that A'2,3 is a circle with centre {Saj — a\)/{5 — 1), where 
5 = |j8i/j82p. If 5 = 1 then K2^}, is a straight line. The three circles are part of a coaxial system. 

Once one has determined the circles Kij the regions Ui can be identified by use of the fol- 
lowing rules. Each component ofU = C\(^i.2 U /Ti 3 U/r2,3) is contained in a single set Ui and 
the unbounded components are contained in ?7i. If C is some part of Ki j and Ui is on one side 
of C then by applying the maximum principle to /,(z) / fj{z) one sees that Uj is on the other side 
of C. If, however, Uk is on one side of C where k ^ i and j then Uk is also on the other side 
of C and the asymptotic spectrum E does not contain any points in C. 

By applying the above rules one can progressively determine which Ut contains each com- 
ponent of U and also eliminate certain curves C C Kij as possible parts of E. Figure [5] portrays 
the eigenvalues of the matrix A with « = 30, a = [2, — 1], j8 = [2,3] and 7=5. The circles 
involved are |z — 2| =2, |z + 1 1 =3 and |z — 4.4| = 3.6. The arcs removed from each circle are 
in accordance with the application of the above rules. The eigenvalue close to 5 is given more 
accurately by A ~ 5.0000104. 
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Figure 5: Eigenvalues of a matrix with h = l and k=l. 

If the two channels have different lengths, say c\n and C2?i, where « — > 00, then only small 
changes are needed. The anti-Stokes curves Ki^ and K\ j, are unchanged but ^2,3 becomes the 
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set of z such that 



z-ay 




z-a2 




Pi 




Pi 





which is only a circle if ci = C2. The general rules for locating the eigenvalues still apply. 

Example 18 (h = 3, k = I). The analysis in the previous example carries through for larger 
values of h. In particular for /j = 3 the reduced matrix has 7 free parameters apart from n and is 



( {z-aiY 



\ 



{z-a2) 
-Pi 



[z-a^] 

-P3 



-P'l' \ 
-P'i 
-P" 

z-r / 



The characteristic polynomial is 

Piz) 



= (z - 7) (z - «i Tiz - «2)"(z - «3)" - Pr iz - a2)"(z - «3)" 
-P^^\z - ai Yiz - as)" - P"^\z - a, )"(z - «2)". 

Assuming that all the a, are distinct, there are 6 anti-Stokes lines, namely 



K\,2 


= {z 


\z- 


«1 


= \Pl\} 




Kl,3 


= {z 


\z- 




= \Pl\} 




Kl,4 


= {z 


\z- 




= m} 




K2,3 


= {z 


\P1\ 


k~ 


«i| = liSil 


k-«2|} 


K2,4 


= {z 


m 




«i| = liSil 


k-«3|} 


K3,4 


= {z 


\P1\ 


\z- 


«3| = \P3\ 


k-«2|} 



These are all circles, or straight lines in degenerate cases. 

In Figure|6l we put n = 30, a = [1,/,—/], j8 = [1,3/2,3/2] and 7=3. One again the eigen- 
values are close to the parts of the anti-Stokes lines obtained by using the rules given above. 



Example 19 {h = k = 2). We consider a graph with two junctions, each containing one point, 
and two channels, each containing n points. Each channel goes from one junction to the other, 
so that the graph has exactly one circuit. The matrix A acts on vectors of length 2n + 2. The 
indices 2n+ l,2n + 2 label the junctions, {l,...,n} label the first channel and {n + I,. . . ,2n} 
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Figure 6: Eigenvalues of a matrix with h = 3 and k=l. 



label the second channel. The matrix A is defined by 

ai ■•• 

"1 



v 




7i 



72 J 



The characteristic polynomial for general n is 
The spectrum converges to 

{z:\{z-a^){z-a2)\ = \m\}, 



(9) 



as n — >^ oo, apart from the possibility of isolated eigenvalues converging to 71 or This curve 
has two points in common with each of the circles |z — £Z,| = |j3;|, provided these circles intersect. 
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We carried out computations for the case n = 30, a = [—1.2, 1.2], j8 = [1.3, 1.3] and 7 = 
[—2,2]. Figure |7] plots the eigenvalues as crosses and the circles \z — 0C;| = |j8,| as dotted curves. 
The curve ^ is simple but non-convex and crosses both circles at ±0.5/. 

The other curves in the figure are the pseudospectral contours of A for £ = 10^'", where 
m = 0, ...,6, the outermost one, for £ = 1, being only partly visible; see |8, 1 1 1 for the definition. 
One deduces from these contours that the eigenvalues are highly unstable; indeed for « = 50 they 
are not easily computable because of rounding errors. Note that the pseudospectral contours are 
related much more strongly to the pair of circles than they are to the eigenvalues. 

If one links the two junctions together weakly by putting 

A{2n+l,2n + 2) = A{2n + 2,2n+ I) = 10"^ 

then the spectrum of A changes radically and approximates the union of the two circles more 
closely as n increases. The case « = 100 is plotted in Figure H] 




-3-2-10123 



Figure 7: Eigenvalues of a matrix with h = 2, k = 2 and n = 30. 
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Figure 8: Eigenvalues of a matrix with h = 2, k = \ and n = 100. 



Example 20(h = 3,k = 2). The {3n + 3) x {3n + 3) matrix 

' i/2 if I <i = j <n, 



-2 


if / = 7 = M + 1 , 


-i/2 


if n + 2 < i = j <2n+ I, 


2 


if i = j = 2n + 2, 


3/2 


if2n + 3<i = j<3n + 2. 





if i = j = 3n + 3, 


1 


if/+l = j, 


1 


if i = 3n + 3, j = 1, 


1 


if i = n+\, j = 3n + 3. 


. 


otherwise. 



has characteristic polynomial 

det(z/ - A) = z{z^- 4)(z - //2)"(z + //2)«(z - 3/2)" -{z-2){z + i/2)\z - 3/2)" - 1 

Its spectrum is depicted in Figure |9] for n = 25. The corresponding graph has 2 junctions and 
3 channels, each of length n. Since h = 3 there are potentially 8 terms in ([T]), but the explicit 
expression of det(z/ — A) shows that only 3 are non-zero. These give rise to the 3 arcs in Figure|9l 
The eigenvalues near ±2 are associated with two of the three diagonal entries at the junctions. 
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-1.5 1 ' ' ' ' ' ' ' ' ' ' 

-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



Figure 9: The spectrum of Example l20l 

6 Localization of Eigenvectors 

In the general model each channel C,- is determined by three parameters, its length e^, the diago- 
nal entries and the non-zero, off-diagonal entries jS^. We may associate with each channel Cr 
the circle 

Sr = {z:\z-a,.\ = \l5r\}. 
If v^"^ is an eigenvector of A^"^ corresponding to the eigenvalue A*^"', it will be shown that the 
coordinates of v^"' corresponding to the channel cl"^ are influenced by the proximity of A^"' to 
the circle Sr. For convenience we label the vertices of channel C^"' by {1, . . . ,ner} so that the 
coordinates of v^"^ corresponding to C,-"' are Vj""*, . . . , vi"|. 

Theorem 21. Let {X^"^} be a sequence such that X^"^^ is an eigenvalue ofA^"\ and let v^"' denote 
any eigenvector corresponding to A^"'. 

1. If there exists 5 > such that 

dist(AW,5,) > 5 (10) 
for all n, then there exists c € (0, 1) such that for each n either 

I ^ i-ll (")l f • 1 

|v- I < c |Vj I for i = 1, . . . ,ner 

or 

|vj,"'^_,-| < c>i"|.| for / = 0, . . . ,ner - 1. 
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2. If there exists C > such that 



dist(AW,5,-) < 



C 



(11) 



for all n, then there exists d > such that for each n either 



,(«) 



or 



(n) 



v|"' =0 for /=!,..., ncr- 

(AW/-AW)vW =0 
- a,)vW = j8,v;ji for / = 1, . . . - 1. 



<d for all /, J G {1, . . . 



Proof The equation 
implies 



If v["^ = then v|"' = for / = 1, . . . So suppose v|"' 7^ 0. Condition (fTOl) implies there 
exists a positive constant c < 1 such that either — a,-)/j3r| < c or |j8r/(A("' — a^)! < c for 
each n. If - ar)/i8r| < c then 



|v- I < c |vj I for i = 1, . . . 



andif |j8,/(AW-a,)| <cthen 



,(«) 



.|_,-| < c'|vl"|| for i = 0,...,ner-\. 

and this proves 1). To prove 2) we observe that condition ([TTI ) implies there exists a > such 
that 



n 



< 1 + 



and so 



for / = 1 , . . . , ncr. Therefore 



n. 



(«)i ^ I («)i ^ /"i , ^ 

V < V; < IH 

\ n 



(i-gr 
(i+f)" 



< 



,(«) 



< 



(1+f)" 
(1-f)" 



for all /,7 G {1, . . . ,ner}. Since (1 it ^ e as « ^ 00 the result follows. 



□ 



If the eigenvalues A'-"' are a positive distance away from all the circles Si then the corre- 
sponding eigenvectors v'") decrease exponentially along all the channels. This suggests that in 
this case the eigenvectors will be concentrated around the junctions. We make this precise as 
follows. 
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Given n and we define the subset C„ jv of tlie grapii (S^") , — >) by 

h 

Cn,N = U ^ ■N<i< ne,- - N}. 

r=l 

We say that a sequence of normalized eigenvectors v^"' are localized around the junctions if for 
all £ > there exists A'^, depending on e, such that ||v(") |c„;v lb < ^ for Note that if 

nmax{ei : I <i<h} <2N 

then C„.A? = 0, so the bound is automatic. The locaUzation condition therefore refers to the 
asymptotic behaviour of A^"^ as « — > oo. 

Theorem 22. Let {v^")} be a sequence of normalized eigenvectors ofA^"^^ with corresponding 
eigenvalues If there exists d > such that for 

dist(A ^"\Sr)'> 5 for all n and all r= l,...,h 

then the v^"^ are localized around the junctions. 

Proof. Let e > and choose A'^ large enough so that 

1-C2 ^ ~h' 

By Theorem |2T] there exists a positive constant c < 1 such that for r = 1 , . . . , /i and all n either 

|v|"''|<c'"' for i = \,...,ner 

or 

|v2-J <c'-i for / = 0,...,«^',-l. 

Therefore 

h ne,—N ^ ^ h ne,—N ^2(^-1) 

r=l i=N r=l i=N ^ ^ 

for all n, as required. □ 
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